{smcl}
{com}{sf}{ul off}{txt}{.-}
      name:  {res}<unnamed>
       {txt}log:  {res}/Users/xpmiki/Dropbox/esaiasson/replicationData/popexp.smcl
  {txt}log type:  {res}smcl
 {txt}opened on:  {res} 9 Sep 2016, 14:42:45
{txt}
{com}. clear
{txt}
{com}. use "populationExperiment.dta"
{txt}
{com}. 
. * randomization check
. reg t i.gender , robust

{txt}Linear regression                               Number of obs     = {res}     1,001
                                                {txt}F(1, 999)         =  {res}     0.98
                                                {txt}Prob > F          = {res}    0.3233
                                                {txt}R-squared         = {res}    0.0010
                                                {txt}Root MSE          =    {res} .50015

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}    Robust
{col 1}           t{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 6}gender {c |}
{space 5}Kvinna  {c |}{col 14}{res}{space 2} .0312468{col 26}{space 2} .0316183{col 37}{space 1}    0.99{col 46}{space 3}0.323{col 54}{space 4}-.0307992{col 67}{space 3} .0932927
{txt}{space 7}_cons {c |}{col 14}{res}{space 2} .4736842{col 26}{space 2} .0224874{col 37}{space 1}   21.06{col 46}{space 3}0.000{col 54}{space 4} .4295563{col 67}{space 3} .5178121
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}{txt}
{com}. reg t i.edu  , robust

{txt}Linear regression                               Number of obs     = {res}     1,001
                                                {txt}F(2, 998)         =  {res}     1.58
                                                {txt}Prob > F          = {res}    0.2072
                                                {txt}R-squared         = {res}    0.0031
                                                {txt}Root MSE          =    {res} .49986

{txt}{hline 23}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 24}{c |}{col 36}    Robust
{col 1}                     t{col 24}{c |}      Coef.{col 36}   Std. Err.{col 48}      t{col 56}   P>|t|{col 64}     [95% Con{col 77}f. Interval]
{hline 23}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 19}edu {c |}
Gymnasium/Motsvarande  {c |}{col 24}{res}{space 2} -.097119{col 36}{space 2} .0687654{col 47}{space 1}   -1.41{col 56}{space 3}0.158{col 64}{space 4}-.2320604{col 77}{space 3} .0378224
{txt}{space 1}Universitet/H�gskola  {c |}{col 24}{res}{space 2}-.0520151{col 36}{space 2} .0686555{col 47}{space 1}   -0.76{col 56}{space 3}0.449{col 64}{space 4}-.1867408{col 77}{space 3} .0827106
{txt}{space 22} {c |}
{space 17}_cons {c |}{col 24}{res}{space 2}  .559322{col 36}{space 2} .0647318{col 47}{space 1}    8.64{col 56}{space 3}0.000{col 64}{space 4} .4322961{col 77}{space 3}  .686348
{txt}{hline 23}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}{txt}
{com}. reg t i.civil  , robust

{txt}Linear regression                               Number of obs     = {res}     1,000
                                                {txt}F(2, 997)         =  {res}     1.80
                                                {txt}Prob > F          = {res}    0.1661
                                                {txt}R-squared         = {res}    0.0036
                                                {txt}Root MSE          =    {res} .49973

{txt}{hline 14}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 15}{c |}{col 27}    Robust
{col 1}            t{col 15}{c |}      Coef.{col 27}   Std. Err.{col 39}      t{col 47}   P>|t|{col 55}     [95% Con{col 68}f. Interval]
{hline 14}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 8}civil {c |}
Sammanboende  {c |}{col 15}{res}{space 2}  .000185{col 27}{space 2} .0402897{col 38}{space 1}    0.00{col 47}{space 3}0.996{col 55}{space 4}-.0788773{col 68}{space 3} .0792473
{txt}{space 7}Annat  {c |}{col 15}{res}{space 2}-.0651773{col 27}{space 2} .0368574{col 38}{space 1}   -1.77{col 47}{space 3}0.077{col 55}{space 4}-.1375042{col 68}{space 3} .0071496
{txt}{space 13} {c |}
{space 8}_cons {c |}{col 15}{res}{space 2} .5085106{col 27}{space 2} .0230946{col 38}{space 1}   22.02{col 47}{space 3}0.000{col 55}{space 4}  .463191{col 68}{space 3} .5538303
{txt}{hline 14}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}{txt}
{com}. reg t i.age , robust

{txt}Linear regression                               Number of obs     = {res}     1,001
                                                {txt}F(3, 997)         =  {res}     0.27
                                                {txt}Prob > F          = {res}    0.8503
                                                {txt}R-squared         = {res}    0.0008
                                                {txt}Root MSE          =    {res} .50069

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}    Robust
{col 1}           t{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}age {c |}
{space 3}30-49 �r  {c |}{col 14}{res}{space 2} .0154314{col 26}{space 2} .0449277{col 37}{space 1}    0.34{col 46}{space 3}0.731{col 54}{space 4}-.0727324{col 67}{space 3} .1035952
{txt}{space 3}50-64 �r  {c |}{col 14}{res}{space 2} .0297315{col 26}{space 2} .0492519{col 37}{space 1}    0.60{col 46}{space 3}0.546{col 54}{space 4}-.0669177{col 67}{space 3} .1263807
{txt}{space 3}65-79 �r  {c |}{col 14}{res}{space 2} .0410058{col 26}{space 2} .0497014{col 37}{space 1}    0.83{col 46}{space 3}0.410{col 54}{space 4}-.0565255{col 67}{space 3} .1385371
{txt}{space 12} {c |}
{space 7}_cons {c |}{col 14}{res}{space 2} .4680851{col 26}{space 2} .0364648{col 37}{space 1}   12.84{col 46}{space 3}0.000{col 54}{space 4} .3965285{col 67}{space 3} .5396417
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}{txt}
{com}. 
. 
. 
. set seed 358395
{txt}
{com}. sem (d1 <-  t  wl p1)(p1 <-  t  wl)   
{res}{txt}
Endogenous variables

{p 0 11 2}Observed:{space 2}{res}d1 p1{p_end}
{txt}
Exogenous variables

{p 0 11 2}Observed:{space 2}{res}t wl{p_end}
{txt}
Fitting target model:

Iteration 0:{space 3}log likelihood = {res:-1649.1573}  
Iteration 1:{space 3}log likelihood = {res:-1649.1573}  (backed up)

Structural equation model{col 49}Number of obs{col 67}= {res}     1,001
{col 1}{txt}Estimation method{col 20}= {res}ml
{col 1}{txt}Log likelihood{col 20}= {res}-1649.1573

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}      OIM
{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}{res}Structural  {col 14}{txt}{c |}
{space 2}{col 3}d1 <-     {col 14}{c |}
{space 10}p1 {c |}{col 14}{res}{space 2}  .467021{col 26}{space 2} .0260961{col 37}{space 1}   17.90{col 46}{space 3}0.000{col 54}{space 4} .4158737{col 67}{space 3} .5181684
{txt}{space 11}t {c |}{col 14}{res}{space 2}-.0335893{col 26}{space 2} .0157188{col 37}{space 1}   -2.14{col 46}{space 3}0.033{col 54}{space 4}-.0643977{col 67}{space 3} -.002781
{txt}{space 10}wl {c |}{col 14}{res}{space 2} .2954601{col 26}{space 2} .0168045{col 37}{space 1}   17.58{col 46}{space 3}0.000{col 54}{space 4} .2625238{col 67}{space 3} .3283964
{txt}{space 7}_cons {c |}{col 14}{res}{space 2} .2466441{col 26}{space 2} .0165188{col 37}{space 1}   14.93{col 46}{space 3}0.000{col 54}{space 4} .2142679{col 67}{space 3} .2790204
{space 2}{txt}{hline 11}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 2}{col 3}p1 <-     {col 14}{c |}
{space 11}t {c |}{col 14}{res}{space 2} .1225325{col 26}{space 2} .0186402{col 37}{space 1}    6.57{col 46}{space 3}0.000{col 54}{space 4} .0859984{col 67}{space 3} .1590666
{txt}{space 10}wl {c |}{col 14}{res}{space 2} .2583145{col 26}{space 2} .0186439{col 37}{space 1}   13.86{col 46}{space 3}0.000{col 54}{space 4} .2217731{col 67}{space 3} .2948559
{txt}{space 7}_cons {c |}{col 14}{res}{space 2} .3880391{col 26}{space 2} .0158071{col 37}{space 1}   24.55{col 46}{space 3}0.000{col 54}{space 4} .3570578{col 67}{space 3} .4190204
{txt}{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
    var(e.d1){c |}{col 14}{res}{space 2} .0589766{col 26}{space 2} .0026362{col 54}{space 4} .0540296{col 67}{space 3} .0643766
{txt}    var(e.p1){c |}{col 14}{res}{space 2} .0865158{col 26}{space 2} .0038672{col 54}{space 4} .0792588{col 67}{space 3} .0944373
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
LR test of model vs. saturated: chi2({res:0})   = {res:     0.00}, Prob > chi2 = {res}     .
{txt}
{com}. estat ic

Akaike's information criterion and Bayesian information criterion

{txt}{hline 13}{c TT}{hline 63}
       Model {c |}        Obs  ll(null)  ll(model)      df         AIC        BIC
{hline 13}{c +}{hline 63}
{ralign 12:.}{col 14}{c |}{res}{col 16}     1,001{col 27}        .{col 38}-1649.157{col 49}     9{col 58} 3316.315{col 69} 3360.493
{txt}{hline 13}{c BT}{hline 63}
{p 15 21 2}
Note: N=Obs used in calculating BIC; see {helpb bic_note:[R] BIC note}.
{p_end}

{com}. estat teffects 
{res}

{txt}Direct effects
{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}      OIM
{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}{res}Structural  {col 14}{txt}{c |}
{space 2}{col 3}d1 <-     {col 14}{c |}
{space 10}p1 {c |}{col 14}{res}{space 2}  .467021{col 26}{space 2} .0260961{col 37}{space 1}   17.90{col 46}{space 3}0.000{col 54}{space 4} .4158737{col 67}{space 3} .5181684
{txt}{space 11}t {c |}{col 14}{res}{space 2}-.0335893{col 26}{space 2} .0157188{col 37}{space 1}   -2.14{col 46}{space 3}0.033{col 54}{space 4}-.0643977{col 67}{space 3} -.002781
{txt}{space 10}wl {c |}{col 14}{res}{space 2} .2954601{col 26}{space 2} .0168045{col 37}{space 1}   17.58{col 46}{space 3}0.000{col 54}{space 4} .2625238{col 67}{space 3} .3283964
{space 2}{txt}{hline 11}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 2}{col 3}p1 <-     {col 14}{c |}
{space 11}t {c |}{col 14}{res}{space 2} .1225325{col 26}{space 2} .0186402{col 37}{space 1}    6.57{col 46}{space 3}0.000{col 54}{space 4} .0859984{col 67}{space 3} .1590666
{txt}{space 10}wl {c |}{col 14}{res}{space 2} .2583145{col 26}{space 2} .0186439{col 37}{space 1}   13.86{col 46}{space 3}0.000{col 54}{space 4} .2217731{col 67}{space 3} .2948559
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}


Indirect effects
{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}      OIM
{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}{res}Structural  {col 14}{txt}{c |}
{space 2}{col 3}d1 <-     {col 14}{c |}
{space 10}p1 {c |}{col 14}{res}{space 2}        0{col 26}{txt}  (no path)
{space 11}t {c |}{col 14}{res}{space 2} .0572253{col 26}{space 2} .0092741{col 37}{space 1}    6.17{col 46}{space 3}0.000{col 54}{space 4} .0390484{col 67}{space 3} .0754021
{txt}{space 10}wl {c |}{col 14}{res}{space 2} .1206383{col 26}{space 2} .0110116{col 37}{space 1}   10.96{col 46}{space 3}0.000{col 54}{space 4}  .099056{col 67}{space 3} .1422206
{space 2}{txt}{hline 11}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 2}{col 3}p1 <-     {col 14}{c |}
{space 11}t {c |}{col 14}{res}{space 2}        0{col 26}{txt}  (no path)
{space 10}wl {c |}{col 14}{res}{space 2}        0{col 26}{txt}  (no path)
{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}


Total effects
{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}      OIM
{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}{res}Structural  {col 14}{txt}{c |}
{space 2}{col 3}d1 <-     {col 14}{c |}
{space 10}p1 {c |}{col 14}{res}{space 2}  .467021{col 26}{space 2} .0260961{col 37}{space 1}   17.90{col 46}{space 3}0.000{col 54}{space 4} .4158737{col 67}{space 3} .5181684
{txt}{space 11}t {c |}{col 14}{res}{space 2} .0236359{col 26}{space 2} .0176816{col 37}{space 1}    1.34{col 46}{space 3}0.181{col 54}{space 4}-.0110194{col 67}{space 3} .0582913
{txt}{space 10}wl {c |}{col 14}{res}{space 2} .4160984{col 26}{space 2} .0176852{col 37}{space 1}   23.53{col 46}{space 3}0.000{col 54}{space 4} .3814361{col 67}{space 3} .4507607
{space 2}{txt}{hline 11}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 2}{col 3}p1 <-     {col 14}{c |}
{space 11}t {c |}{col 14}{res}{space 2} .1225325{col 26}{space 2} .0186402{col 37}{space 1}    6.57{col 46}{space 3}0.000{col 54}{space 4} .0859984{col 67}{space 3} .1590666
{txt}{space 10}wl {c |}{col 14}{res}{space 2} .2583145{col 26}{space 2} .0186439{col 37}{space 1}   13.86{col 46}{space 3}0.000{col 54}{space 4} .2217731{col 67}{space 3} .2948559
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. 
. *** indirect and direct - winning/losing
. cap program drop indireff
{txt}
{com}. program indireff, rclass
{txt}  1{com}. sem (d1 <-  t  wl p1)(p1 <-  t  wl)  
{txt}  2{com}.   estat teffects
{txt}  3{com}.   mat bi = r(indirect)
{txt}  4{com}.   mat bt = r(total)
{txt}  5{com}.   return scalar indir  = el(bi,1,3)
{txt}  6{com}.   return scalar total  = el(bt,1,3)
{txt}  7{com}. end 
{txt}
{com}. 
. set seed 358395
{txt}
{com}. bootstrap r(indir) r(total), reps(1000): indireff 
{txt}(running indireff on estimation sample)

Bootstrap replications ({res}1000{txt})
{hline 4}{c +}{hline 3} 1 {hline 3}{c +}{hline 3} 2 {hline 3}{c +}{hline 3} 3 {hline 3}{c +}{hline 3} 4 {hline 3}{c +}{hline 3} 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000
{res}
{txt}Bootstrap results{col 49}Number of obs{col 67}= {res}     1,001
{txt}{col 49}Replications{col 67}= {res}     1,000
{p2colset 7 17 21 2}{...}

{txt}{p2col :command:}indireff{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_1:}{res:r(indir)}{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_2:}{res:r(total)}{p_end}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}   Observed{col 26}   Bootstrap{col 54}         Norm{col 67}al-based
{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 7}_bs_1 {c |}{col 14}{res}{space 2} .1206383{col 26}{space 2} .0121495{col 37}{space 1}    9.93{col 46}{space 3}0.000{col 54}{space 4} .0968258{col 67}{space 3} .1444508
{txt}{space 7}_bs_2 {c |}{col 14}{res}{space 2} .4160984{col 26}{space 2} .0178393{col 37}{space 1}   23.32{col 46}{space 3}0.000{col 54}{space 4}  .381134{col 67}{space 3} .4510628
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. estat bootstrap, percentile bc  

{txt}Bootstrap results{col 49}Number of obs{col 67}= {res}     1,001
{txt}{col 49}Replications{col 67}= {res}      1000
{p2colset 7 17 21 2}{...}

{txt}{p2col :command:}indireff{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_1:}{res:r(indir)}{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_2:}{res:r(total)}{p_end}

{col 1}{text}{hline 13}{c TT}{hline 12}{hline 11}{hline 13}{hline 12}{hline 10}{hline 6}
{col 14}{text}{c |}    Observed{col 38}    Bootstrap
{col 14}{text}{c |}       Coef.{col 27}       Bias{col 38}    Std. Err.{col 51}  [95% Conf. Interval]
{col 1}{text}{hline 13}{c +}{hline 12}{hline 11}{hline 13}{hline 12}{hline 10}{hline 6}
{col 1}{text}       _bs_1{col 14}{c |}{result}{space 2} .12063831{col 27}{space 2} .0004255{col 38}{space 2} .01214947{col 51}{space 2} .0987132{col 63}{space 1} .1459039{col 73}{text}   (P)
{col 14}{text}{c |}{col 51}{result}{space 2} .0984599{col 63}{space 1} .1453353{col 73}{text}  (BC)
{col 1}{text}       _bs_2{col 14}{c |}{result}{space 2} .41609841{col 27}{space 2}-.0003869{col 38}{space 2} .01783932{col 51}{space 2} .3794943{col 63}{space 1} .4509158{col 73}{text}   (P)
{col 14}{text}{c |}{col 51}{result}{space 2} .3790663{col 63}{space 1} .4507539{col 73}{text}  (BC)
{col 1}{text}{hline 13}{c BT}{hline 12}{hline 11}{hline 13}{hline 12}{hline 10}{hline 6}
{col 0}(P){col 8}percentile confidence interval
{col 0}(BC){col 8}bias-corrected confidence interval

{com}. 
. *** indirect and direct - treatment
. cap program drop indireff
{txt}
{com}. program indireff, rclass
{txt}  1{com}. sem (d1 <- wl  t  p1)(p1 <- wl  t )  
{txt}  2{com}.   estat teffects
{txt}  3{com}.   mat bi = r(indirect)
{txt}  4{com}.   mat bt = r(total)
{txt}  5{com}.   return scalar indirect  = el(bi,1,3)
{txt}  6{com}.   return scalar total  = el(bt,1,3)
{txt}  7{com}. end
{txt}
{com}. 
. set seed 358395
{txt}
{com}. bootstrap r(indirect) r(total), reps(1000): indireff
{txt}(running indireff on estimation sample)

Bootstrap replications ({res}1000{txt})
{hline 4}{c +}{hline 3} 1 {hline 3}{c +}{hline 3} 2 {hline 3}{c +}{hline 3} 3 {hline 3}{c +}{hline 3} 4 {hline 3}{c +}{hline 3} 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000
{res}
{txt}Bootstrap results{col 49}Number of obs{col 67}= {res}     1,001
{txt}{col 49}Replications{col 67}= {res}     1,000
{p2colset 7 17 21 2}{...}

{txt}{p2col :command:}indireff{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_1:}{res:r(indirect)}{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_2:}{res:r(total)}{p_end}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}   Observed{col 26}   Bootstrap{col 54}         Norm{col 67}al-based
{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 7}_bs_1 {c |}{col 14}{res}{space 2} .0572253{col 26}{space 2} .0100498{col 37}{space 1}    5.69{col 46}{space 3}0.000{col 54}{space 4} .0375279{col 67}{space 3} .0769226
{txt}{space 7}_bs_2 {c |}{col 14}{res}{space 2} .0236359{col 26}{space 2} .0175541{col 37}{space 1}    1.35{col 46}{space 3}0.178{col 54}{space 4}-.0107695{col 67}{space 3} .0580413
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. estat bootstrap, percentile bc  

{txt}Bootstrap results{col 49}Number of obs{col 67}= {res}     1,001
{txt}{col 49}Replications{col 67}= {res}      1000
{p2colset 7 17 21 2}{...}

{txt}{p2col :command:}indireff{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_1:}{res:r(indirect)}{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_2:}{res:r(total)}{p_end}

{col 1}{text}{hline 13}{c TT}{hline 12}{hline 11}{hline 13}{hline 12}{hline 10}{hline 6}
{col 14}{text}{c |}    Observed{col 38}    Bootstrap
{col 14}{text}{c |}       Coef.{col 27}       Bias{col 38}    Std. Err.{col 51}  [95% Conf. Interval]
{col 1}{text}{hline 13}{c +}{hline 12}{hline 11}{hline 13}{hline 12}{hline 10}{hline 6}
{col 1}{text}       _bs_1{col 14}{c |}{result}{space 2} .05722525{col 27}{space 2} .0002342{col 38}{space 2} .01004984{col 51}{space 2} .0383927{col 63}{space 1} .0782295{col 73}{text}   (P)
{col 14}{text}{c |}{col 51}{result}{space 2} .0381328{col 63}{space 1} .0770452{col 73}{text}  (BC)
{col 1}{text}       _bs_2{col 14}{c |}{result}{space 2} .02363592{col 27}{space 2} .0002156{col 38}{space 2} .01755409{col 51}{space 2}-.0103882{col 63}{space 1} .0587565{col 73}{text}   (P)
{col 14}{text}{c |}{col 51}{result}{space 2}-.0097948{col 63}{space 1} .0598807{col 73}{text}  (BC)
{col 1}{text}{hline 13}{c BT}{hline 12}{hline 11}{hline 13}{hline 12}{hline 10}{hline 6}
{col 0}(P){col 8}percentile confidence interval
{col 0}(BC){col 8}bias-corrected confidence interval

{com}. 
. 
. ********
. ******** ROBUSTNESS CHECKS
. ********
. 
. set seed 358395
{txt}
{com}. sem (d2 <-  t  wl p2)(p2 <-  t  wl)   
{res}{txt}
Endogenous variables

{p 0 11 2}Observed:{space 2}{res}d2 p2{p_end}
{txt}
Exogenous variables

{p 0 11 2}Observed:{space 2}{res}t wl{p_end}
{txt}
Fitting target model:

Iteration 0:{space 3}log likelihood = {res:-1661.8417}  
Iteration 1:{space 3}log likelihood = {res:-1661.8417}  

Structural equation model{col 49}Number of obs{col 67}= {res}     1,001
{col 1}{txt}Estimation method{col 20}= {res}ml
{col 1}{txt}Log likelihood{col 20}= {res}-1661.8417

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}      OIM
{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}{res}Structural  {col 14}{txt}{c |}
{space 2}{col 3}d2 <-     {col 14}{c |}
{space 10}p2 {c |}{col 14}{res}{space 2} .3844835{col 26}{space 2} .0270252{col 37}{space 1}   14.23{col 46}{space 3}0.000{col 54}{space 4} .3315151{col 67}{space 3} .4374519
{txt}{space 11}t {c |}{col 14}{res}{space 2}-.0303986{col 26}{space 2} .0159404{col 37}{space 1}   -1.91{col 46}{space 3}0.057{col 54}{space 4}-.0616412{col 67}{space 3} .0008441
{txt}{space 10}wl {c |}{col 14}{res}{space 2} .2949702{col 26}{space 2} .0167793{col 37}{space 1}   17.58{col 46}{space 3}0.000{col 54}{space 4} .2620834{col 67}{space 3}  .327857
{txt}{space 7}_cons {c |}{col 14}{res}{space 2} .3192513{col 26}{space 2} .0170838{col 37}{space 1}   18.69{col 46}{space 3}0.000{col 54}{space 4} .2857677{col 67}{space 3} .3527349
{space 2}{txt}{hline 11}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 2}{col 3}p2 <-     {col 14}{c |}
{space 11}t {c |}{col 14}{res}{space 2} .0881843{col 26}{space 2} .0184334{col 37}{space 1}    4.78{col 46}{space 3}0.000{col 54}{space 4} .0520555{col 67}{space 3} .1243131
{txt}{space 10}wl {c |}{col 14}{res}{space 2} .2126527{col 26}{space 2} .0184371{col 37}{space 1}   11.53{col 46}{space 3}0.000{col 54}{space 4} .1765167{col 67}{space 3} .2487887
{txt}{space 7}_cons {c |}{col 14}{res}{space 2} .3937132{col 26}{space 2} .0156317{col 37}{space 1}   25.19{col 46}{space 3}0.000{col 54}{space 4} .3630756{col 67}{space 3} .4243508
{txt}{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
    var(e.d2){c |}{col 14}{res}{space 2} .0618553{col 26}{space 2} .0027649{col 54}{space 4} .0566668{col 67}{space 3} .0675188
{txt}    var(e.p2){c |}{col 14}{res}{space 2} .0846068{col 26}{space 2} .0037818{col 54}{space 4} .0775099{col 67}{space 3} .0923534
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
LR test of model vs. saturated: chi2({res:0})   = {res:     0.00}, Prob > chi2 = {res}     .
{txt}
{com}. outreg2 using popExp2.doc, replace ctitle(1) dec(3)
{txt}{stata `"shellout using `"popExp2.doc"'"':popExp2.doc}
{browse `"/Users/xpmiki/Dropbox/esaiasson/replicationData"' :dir}{com} : {txt}{stata `"seeout using "popExp2.txt""':seeout}

{com}. estat ic

Akaike's information criterion and Bayesian information criterion

{txt}{hline 13}{c TT}{hline 63}
       Model {c |}        Obs  ll(null)  ll(model)      df         AIC        BIC
{hline 13}{c +}{hline 63}
{ralign 12:.}{col 14}{c |}{res}{col 16}     1,001{col 27}        .{col 38}-1661.842{col 49}     9{col 58} 3341.683{col 69} 3385.862
{txt}{hline 13}{c BT}{hline 63}
{p 15 21 2}
Note: N=Obs used in calculating BIC; see {helpb bic_note:[R] BIC note}.
{p_end}

{com}. estat teffects 
{res}

{txt}Direct effects
{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}      OIM
{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}{res}Structural  {col 14}{txt}{c |}
{space 2}{col 3}d2 <-     {col 14}{c |}
{space 10}p2 {c |}{col 14}{res}{space 2} .3844835{col 26}{space 2} .0270252{col 37}{space 1}   14.23{col 46}{space 3}0.000{col 54}{space 4} .3315151{col 67}{space 3} .4374519
{txt}{space 11}t {c |}{col 14}{res}{space 2}-.0303986{col 26}{space 2} .0159404{col 37}{space 1}   -1.91{col 46}{space 3}0.057{col 54}{space 4}-.0616412{col 67}{space 3} .0008441
{txt}{space 10}wl {c |}{col 14}{res}{space 2} .2949702{col 26}{space 2} .0167793{col 37}{space 1}   17.58{col 46}{space 3}0.000{col 54}{space 4} .2620834{col 67}{space 3}  .327857
{space 2}{txt}{hline 11}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 2}{col 3}p2 <-     {col 14}{c |}
{space 11}t {c |}{col 14}{res}{space 2} .0881843{col 26}{space 2} .0184334{col 37}{space 1}    4.78{col 46}{space 3}0.000{col 54}{space 4} .0520555{col 67}{space 3} .1243131
{txt}{space 10}wl {c |}{col 14}{res}{space 2} .2126527{col 26}{space 2} .0184371{col 37}{space 1}   11.53{col 46}{space 3}0.000{col 54}{space 4} .1765167{col 67}{space 3} .2487887
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}


Indirect effects
{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}      OIM
{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}{res}Structural  {col 14}{txt}{c |}
{space 2}{col 3}d2 <-     {col 14}{c |}
{space 10}p2 {c |}{col 14}{res}{space 2}        0{col 26}{txt}  (no path)
{space 11}t {c |}{col 14}{res}{space 2} .0339054{col 26}{space 2} .0074773{col 37}{space 1}    4.53{col 46}{space 3}0.000{col 54}{space 4} .0192502{col 67}{space 3} .0485606
{txt}{space 10}wl {c |}{col 14}{res}{space 2} .0817615{col 26}{space 2} .0091257{col 37}{space 1}    8.96{col 46}{space 3}0.000{col 54}{space 4} .0638754{col 67}{space 3} .0996475
{space 2}{txt}{hline 11}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 2}{col 3}p2 <-     {col 14}{c |}
{space 11}t {c |}{col 14}{res}{space 2}        0{col 26}{txt}  (no path)
{space 10}wl {c |}{col 14}{res}{space 2}        0{col 26}{txt}  (no path)
{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}


Total effects
{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}      OIM
{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}{res}Structural  {col 14}{txt}{c |}
{space 2}{col 3}d2 <-     {col 14}{c |}
{space 10}p2 {c |}{col 14}{res}{space 2} .3844835{col 26}{space 2} .0270252{col 37}{space 1}   14.23{col 46}{space 3}0.000{col 54}{space 4} .3315151{col 67}{space 3} .4374519
{txt}{space 11}t {c |}{col 14}{res}{space 2} .0035068{col 26}{space 2} .0172814{col 37}{space 1}    0.20{col 46}{space 3}0.839{col 54}{space 4}-.0303641{col 67}{space 3} .0373778
{txt}{space 10}wl {c |}{col 14}{res}{space 2} .3767317{col 26}{space 2} .0172849{col 37}{space 1}   21.80{col 46}{space 3}0.000{col 54}{space 4} .3428539{col 67}{space 3} .4106094
{space 2}{txt}{hline 11}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 2}{col 3}p2 <-     {col 14}{c |}
{space 11}t {c |}{col 14}{res}{space 2} .0881843{col 26}{space 2} .0184334{col 37}{space 1}    4.78{col 46}{space 3}0.000{col 54}{space 4} .0520555{col 67}{space 3} .1243131
{txt}{space 10}wl {c |}{col 14}{res}{space 2} .2126527{col 26}{space 2} .0184371{col 37}{space 1}   11.53{col 46}{space 3}0.000{col 54}{space 4} .1765167{col 67}{space 3} .2487887
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. 
. *** indirect and direct - winning/losing
. cap program drop indireff
{txt}
{com}. program indireff, rclass
{txt}  1{com}. sem (d2 <-  t  wl p2)(p2 <-  t  wl)  
{txt}  2{com}.   estat teffects
{txt}  3{com}.   mat bi = r(indirect)
{txt}  4{com}.   mat bt = r(total)
{txt}  5{com}.   return scalar indir  = el(bi,1,3)
{txt}  6{com}.   return scalar total  = el(bt,1,3)
{txt}  7{com}. end 
{txt}
{com}. 
. set seed 358395
{txt}
{com}. bootstrap r(indir) r(total), reps(1000): indireff 
{txt}(running indireff on estimation sample)

Bootstrap replications ({res}1000{txt})
{hline 4}{c +}{hline 3} 1 {hline 3}{c +}{hline 3} 2 {hline 3}{c +}{hline 3} 3 {hline 3}{c +}{hline 3} 4 {hline 3}{c +}{hline 3} 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000
{res}
{txt}Bootstrap results{col 49}Number of obs{col 67}= {res}     1,001
{txt}{col 49}Replications{col 67}= {res}     1,000
{p2colset 7 17 21 2}{...}

{txt}{p2col :command:}indireff{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_1:}{res:r(indir)}{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_2:}{res:r(total)}{p_end}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}   Observed{col 26}   Bootstrap{col 54}         Norm{col 67}al-based
{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 7}_bs_1 {c |}{col 14}{res}{space 2} .0817615{col 26}{space 2} .0099492{col 37}{space 1}    8.22{col 46}{space 3}0.000{col 54}{space 4} .0622614{col 67}{space 3} .1012615
{txt}{space 7}_bs_2 {c |}{col 14}{res}{space 2} .3767317{col 26}{space 2} .0175899{col 37}{space 1}   21.42{col 46}{space 3}0.000{col 54}{space 4} .3422561{col 67}{space 3} .4112072
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. estat bootstrap, percentile bc  

{txt}Bootstrap results{col 49}Number of obs{col 67}= {res}     1,001
{txt}{col 49}Replications{col 67}= {res}      1000
{p2colset 7 17 21 2}{...}

{txt}{p2col :command:}indireff{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_1:}{res:r(indir)}{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_2:}{res:r(total)}{p_end}

{col 1}{text}{hline 13}{c TT}{hline 12}{hline 11}{hline 13}{hline 12}{hline 10}{hline 6}
{col 14}{text}{c |}    Observed{col 38}    Bootstrap
{col 14}{text}{c |}       Coef.{col 27}       Bias{col 38}    Std. Err.{col 51}  [95% Conf. Interval]
{col 1}{text}{hline 13}{c +}{hline 12}{hline 11}{hline 13}{hline 12}{hline 10}{hline 6}
{col 1}{text}       _bs_1{col 14}{c |}{result}{space 2} .08176146{col 27}{space 2}  .000251{col 38}{space 2} .00994919{col 51}{space 2}  .063342{col 63}{space 1} .1023574{col 73}{text}   (P)
{col 14}{text}{c |}{col 51}{result}{space 2} .0632616{col 63}{space 1}  .102248{col 73}{text}  (BC)
{col 1}{text}       _bs_2{col 14}{c |}{result}{space 2} .37673166{col 27}{space 2}-.0007051{col 38}{space 2} .01758991{col 51}{space 2} .3419997{col 63}{space 1} .4114601{col 73}{text}   (P)
{col 14}{text}{c |}{col 51}{result}{space 2}  .343565{col 63}{space 1} .4125254{col 73}{text}  (BC)
{col 1}{text}{hline 13}{c BT}{hline 12}{hline 11}{hline 13}{hline 12}{hline 10}{hline 6}
{col 0}(P){col 8}percentile confidence interval
{col 0}(BC){col 8}bias-corrected confidence interval

{com}. 
. *** indirect and direct - treatment
. cap program drop indireff
{txt}
{com}. program indireff, rclass
{txt}  1{com}. sem (d2 <- wl  t  p2)(p2 <- wl  t )  
{txt}  2{com}.   estat teffects
{txt}  3{com}.   mat bi = r(indirect)
{txt}  4{com}.   mat bt = r(total)
{txt}  5{com}.   return scalar indirect  = el(bi,1,3)
{txt}  6{com}.   return scalar total  = el(bt,1,3)
{txt}  7{com}. end
{txt}
{com}. 
. set seed 358395
{txt}
{com}. bootstrap r(indirect) r(total), reps(1000): indireff
{txt}(running indireff on estimation sample)

Bootstrap replications ({res}1000{txt})
{hline 4}{c +}{hline 3} 1 {hline 3}{c +}{hline 3} 2 {hline 3}{c +}{hline 3} 3 {hline 3}{c +}{hline 3} 4 {hline 3}{c +}{hline 3} 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000
{res}
{txt}Bootstrap results{col 49}Number of obs{col 67}= {res}     1,001
{txt}{col 49}Replications{col 67}= {res}     1,000
{p2colset 7 17 21 2}{...}

{txt}{p2col :command:}indireff{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_1:}{res:r(indirect)}{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_2:}{res:r(total)}{p_end}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}   Observed{col 26}   Bootstrap{col 54}         Norm{col 67}al-based
{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 7}_bs_1 {c |}{col 14}{res}{space 2} .0339054{col 26}{space 2} .0080041{col 37}{space 1}    4.24{col 46}{space 3}0.000{col 54}{space 4} .0182177{col 67}{space 3} .0495931
{txt}{space 7}_bs_2 {c |}{col 14}{res}{space 2} .0035068{col 26}{space 2} .0177007{col 37}{space 1}    0.20{col 46}{space 3}0.843{col 54}{space 4} -.031186{col 67}{space 3} .0381996
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. estat bootstrap, percentile bc  

{txt}Bootstrap results{col 49}Number of obs{col 67}= {res}     1,001
{txt}{col 49}Replications{col 67}= {res}      1000
{p2colset 7 17 21 2}{...}

{txt}{p2col :command:}indireff{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_1:}{res:r(indirect)}{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_2:}{res:r(total)}{p_end}

{col 1}{text}{hline 13}{c TT}{hline 12}{hline 11}{hline 13}{hline 12}{hline 10}{hline 6}
{col 14}{text}{c |}    Observed{col 38}    Bootstrap
{col 14}{text}{c |}       Coef.{col 27}       Bias{col 38}    Std. Err.{col 51}  [95% Conf. Interval]
{col 1}{text}{hline 13}{c +}{hline 12}{hline 11}{hline 13}{hline 12}{hline 10}{hline 6}
{col 1}{text}       _bs_1{col 14}{c |}{result}{space 2}  .0339054{col 27}{space 2} .0002554{col 38}{space 2} .00800409{col 51}{space 2} .0197407{col 63}{space 1} .0510565{col 73}{text}   (P)
{col 14}{text}{c |}{col 51}{result}{space 2} .0196316{col 63}{space 1} .0501077{col 73}{text}  (BC)
{col 1}{text}       _bs_2{col 14}{c |}{result}{space 2} .00350684{col 27}{space 2} .0000824{col 38}{space 2} .01770073{col 51}{space 2}-.0314683{col 63}{space 1} .0370955{col 73}{text}   (P)
{col 14}{text}{c |}{col 51}{result}{space 2}-.0328646{col 63}{space 1} .0352743{col 73}{text}  (BC)
{col 1}{text}{hline 13}{c BT}{hline 12}{hline 11}{hline 13}{hline 12}{hline 10}{hline 6}
{col 0}(P){col 8}percentile confidence interval
{col 0}(BC){col 8}bias-corrected confidence interval

{com}. 
. log close
      {txt}name:  {res}<unnamed>
       {txt}log:  {res}/Users/xpmiki/Dropbox/esaiasson/replicationData/popexp.smcl
  {txt}log type:  {res}smcl
 {txt}closed on:  {res} 9 Sep 2016, 14:46:31
{txt}{.-}
{smcl}
{txt}{sf}{ul off}